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We propose to use the thermal model with the multi-component hard-core radii to describe 
the hadron yield ratios from the low AGS to the highest RHIC energies. It is demonstrated 
that the variation of the hard-core radii of pions and kaons enable us to drastically improve 
the fit quality of the measured mid-rapidity data and for the first time to completely describe 
the Strangeness Horn behavior as the function of the energy of collision without spoiling the 
fit quality of other ratios. The best global fit is found for the vanishing hard-core radius of 
pions and for the hard-core radius of kaons being equal to 0.35 fm, whereas the hard-core 
radius of all other mesons is fixed to 0.3 fm and that one of baryons is fixed to 0.5 fm. 
It is argued that the multi-component hadron resonance gas model opens us a principal 
possibility to determine the second virial coefficients of hadron-hadron interaction. 



I. INTRODUCTION 

The hadron resonance gas model 1 [TJ [2] is the only theoretical tool allowing us to extract information 
about the chemical freeze-out (FO) stage of the relativistic heavy ion collisions. Although its systematic 
application to the experimental data description began about fifteen years ago [3J, many features of this 
model are not well studied [U [5]. Thus, very recently in a critical analysis of the hadron resonance gas 
model [5] it was shown that for the description of the hadron multiplicities the baryon charge conservation 
and the isospin conservation, used in one of the most successful versions of this model pQ, should be 
essentially modified, whereas for the description of the hadron yield ratios these conservation laws are not 
necessary at all. Although the discussion about the reliable chemical FO criterion has a long history [IJ|6], 
only recently it was demonstrated that none of the previously suggested chemical FO criteria, including 
the most popular one of constant energy per particle E/N ~ 1.1 GeV [6], is robust [5], if the realistic 
particle table with the hadron masses up to 2.5 GeV is used. At the same time in [5] it was shown that 
despite an essential difference with the approach used in \1] , the both versions of the hadron resonance gas 
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model demonstrate almost the same value 7.18 for the entropy per particle at chemical FO. Thus, it turns 
out that the criterion of the constant entropy per particle at chemical FO is, indeed, a reliable one. It is 
interesting that the constant entropy per particle at chemical FO was also found in [7], but the way of hard 
core repulsion used in this work is too different from the traditional one used in the hadron resonance gas 
model [HO [5] and, hence, in contrast to the results of [HO [5], the model used in [7] leads to a simultaneous 
fulfillment of a few chemical FO criteria. 

One of the traditional difficulties of the hadron resonance gas model is related to the Strangeness Horn 
description which up to now is far from being satisfactory, although very different formulations of the 
thermal model are used for this purpose [U [5j [8] . Note that the principal importance to improve the 
Strangeness Horn description can be easily understood from the fact that just the non-monotonic behavior 
of the K + /tt + ratio as the function of the center of mass energy of collision is often claimed to be one of a 
few existing signals of the onset on deconfinement [9l411| . The previous attempts [TJ [TJ |8] to describe the 
Strangeness Horn behavior without spoiling the quality of other particle ratios fit and the thorough analysis 
performed in |5j led us to a conclusion that further improvement of the hadron resonance gas model can be 
achieved, if we consider the pion and kaon hard-core radii, as an independent fitting parameters. Evidently, 
this would allow us to have two additional fitting parameters to improve the fit quality of the Strangeness 
Horn without spoiling the other hadron yield ratios. The physical idea behind such an approach is that 
the hadronic hard-core radii are the effective parameters which include the contributions of the repulsion 
and attraction. Since the parameters of hadron-hadron interaction are, generally speaking, individual for 
each kind of hadrons, then each kind of hadrons can have its own hard-core radius. 

The work is organized as follows. Section II contains the main equations of the model. The results are 
discussed in Section III, while the last Section contains our conclusions. 

II. MULTI-COMPONENT HADRON GAS AND HARD-CORE RADII 

The hadron resonance gas model is a successful compromise between the right choice of the physically 
relevant degrees of freedom and the simple parameterization of their interaction. Its theoretical justification 
is based on a simple fact found rather long ago |12j that for temperatures below 170 MeV the interacting 
mixture of stable hadrons and their resonances behaves as the mixture of nearly ideal gases of stable 
particles which in this case include both the hadrons and the resonances taken with their averaged masses. 
The reason for such a behavior is nearly a complete cancellation between the attraction and repulsion 
contributions. The resulting deviation from the ideal gas (a weak repulsion) is usually attributed to the 
second virial coefficients bij defined for the hadrons of z-th and j-th. kinds. Since the equations of state of 
the hadron resonance gas have the Van der Waals type repulsion, the coefficients bij are called as excluded 
volumes. Up to now the hadron resonance gas model employed only two basic parameters related to 
hadron-hadron repulsion: the model of [3] had one excluded volume for pions and another for all other 
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hadrons, whereas the model [13J suggested to consider one common hard-core radius for mesons R m and 
another hard-core radius for baryons Rb- However, none of the models developed in [31 [13] were correct 
[3], since they did not include the second virial coefficient of the crossed type b^ with i ^ j, i.e. between 
the hadrons of different kinds. On the one hand the realization of this fact led to a systematic description 
of the data with a single hard-core radius for all hadrons [I], and on the other hand it also led to the 
development of the multi-component Van der Waals gas models [31 El IE]. 

The success of the hadron resonance gas model, the one pQ or two component [H in the data 
description may look surprising at the first glance, but this is not just a single example of a simple statistical 
model that is able to efficiently account for the complex features of interaction between the constituents. 
One should remember, although the interaction between the clusters of many molecules in the real gases 
or interaction between the nuclear fragments is no less, but more complex than interaction of hadrons, 
the successful statistical models for such systems are well known [16H20] . and, nevertheless, these models 
are able not only to describe the low density states of real gases |16H18j or that ones of nuclear fragments 
|19} I20| . but they are able to successfully model the condensation of these gases into the corresponding 
liquids at rather high densities. And one should also remember that these models employ a few statistical 
parameters only and use rather simple, but physically adequate (!) parameterization for the many-body 
effects. Thus, the whole point is that in all these successful examples |16H19j the employed parameters 
which characterize the interaction are effective from the very beginning and only at very low densities 
the models recover the virial expansion up to the second virial coefficients. Therefore, one should not be 
surprised that the second virial coefficients bij of the hadron resonance gas are some effective parameters 
which account for (to large extent) a cancellation of the attractive and repulsive contributions, and which, 
in principle, could be individual characteristics for each hadronic pair. However, the overall success of the 
hadron resonance gas model evidences that the number of independent parameters should be essentially 
smaller then the number of hadron types. Thus, below we demonstrate that the available data favor just 
the set of the excluded volumes defined as bij = ^-{Ri + Rj) 3 via the hard-core radii of pions R^, kaons 
Rk, baryons R^ and the radius for all other mesons R m . In what follows we give the main equations of 
the multi-component formulation referring to [H [5] for a detailed derivation. 

Consider the Boltzmann gas of N hadron species in a volume V that has the temperature T, the 
baryonic chemical potential fj,g, the strange chemical potential us and the chemical potential of the isospin 
third component The system pressure p and the K-th charge density nf- (K £ {B, S, 13}) of the i-th 
hadron sort are given by the expressions (B denotes a symmetric matrix of the second virial coefficients 
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Here the full chemical potential of the i-th hadron sort = Qf /j,b + Qf Ms + Qf^n is expressed in terms 
of the corresponding charges Qf and their chemical potentials, <j>i(T) denotes the thermal particle density 
of the i-th hadron sort of mass m, and degeneracy gi, and £ T denotes the row of variables £j. 

In a special case when all the elements of the second virial coefficients matrix are equal fry = Vq Eqs. 
Q-Q, evidently, reproduce the one component model with the pressure 
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which defines the particle density of i-th kind of hadron as rij 
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The latter shows 



that the ratios of two particle densities defined by ^ match that ones of the mixture of the corresponding 
ideal gases for an arbitrary value of vq, while the particle densities themselves may essentially differ from 
the particle densities of the ideal gas. 

It is known that the resonance width is important at low temperatures pQ. Similarly to [I], the width 
Ti of the resonance of mean mass rrn is modeled by replacing the Boltzmann distribution function in the 
particle thermal density ([2]) by its average over the Breit-Wigner mass distribution as 
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where Mq is the dominant decay channel mass. Such a substitution provides a simple, but reliable approx- 
imation to account for the resonance width. 

The contribution of the resonance decays is accounted for as usual: the total density of hadron X 



consists of the thermal part and the decay ones: 
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(5) 



where Br{Y — > X) is the decay branching of the Y-th hadron into the hadron X. The masses, the widths and 
the strong decay branchings of all hadrons were taken from the particle tables used by the thermodynamic 
code THERMUS EH- 
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The strange charge conservation completes the list of equations used. Since in strong decays the 
strangeness is conserved, then it is sufficient to impose the vanishing of the total strangeness for ther- 
mal densities at chemical FO, i.e. to determine the strange chemical potential (jls from the equation 

N 

£ nf = 0. As it was shown recently [3] the baryonic charge and isospin conservation laws should not be 

i=l 

imposed to fit the hadron multipllicities since they lead to unphysically huge FO volumes. Therefore, in 
this work for the data at given energy of collision we use the following fitting parameters: temperature 
T, baryonic chemical potential \xb an d the chemical potential of the third projection of isospin Note 
that such a procedure is completely consistent with fitting the hadron multiplicities instead of hadron yield 
ratios [5], the main difference is only that to fit the hadron multiplicities one has to use the chemical FO 
volume V as an additional parameter. As it was explained earlier the global fitting parameters are the 
hard-core radii of pions R n , kaons Rk, baryons Rf, and that one for all other mesons R m . In addition, to 
demonstrate the pure effect of the radii variation we do not include any strangeness suppression factor into 
simulations. Then one should expect some minor problems with the description of multi-strange baryons. 

III. RESULTS 

The recent comprehensive analysis [5] performed for different hard-core radii of all mesons R m and all 
baryons Rf, clearly showed us that the good description of the data can be achieved for many pairs of these 
radii. However, the ratios are more stable during the fitting, if R m = 0.3 fm and Rb = 0.5 fm. Right these 
values of hard-core radii were fixed and then we fitted the data by the x 2 / do f -criterion for different values 
of the pion R^ and kaon Rk hard-core radii taken below 0.5 fm each. The minimal value of X 2 /dof ~ 1.018 
for the energies in the range ^snn = 2.7, 3.3, 3.8, 4.3, 4.9, 6.3, 7.6, 8.8, 12, 17, 130, 200 GeV (for details see 
below) was obtained for R v = fm and Rk = 0.35 fm. 

Since in the present approach there is no principle difference between fitting the absolute hadron yields 
at mid-rapidity or their ratios, we prefer to fit ratios in order to reduce the volume of numerical efforts. 
In our choice of the data sets we basically followed Ref. [IJ. Thus, at the AGS energy range of collisions 
(y/SNN = 2.7 — 4.9 GeV) the data are available for the kinetic beam energies from 2 to 10.7 AGeV. For 
the beam energies 2, 4, 6 and 8 AGeV there are only a few data points available: the yields for pions 
[221 123] . for protons |24|, [25], for kaons |23] (except for 2 AGeV), for A hyperons the integrated over 4tt 
data are available [25]. For the beam energy 6 AGeV there exist the H~ hyperon data integrated over 
Air geometry |27j . However, the data for the A and H~ hyperons have to be corrected [I], and instead of 
the raw experimental data we used their corrected values of Ref. [1] . For the highest AGS center of mass 
energy ^snn = 4.9 GeV (or the beam energy 10.7 AGeV) in addition to the mentioned data for pions, 
(anti)protons and kaons there exist data for <f> meson [28], for A hyperon ||29] and A hyperon |30j . 

As one can see from the left panel of Fig. [I] the quality of the fit achieved for ^/~snn = 4.9 GeV is 
extremely good even for A/A and A/ir~ ratios, i.e. for the most problematic ratios of [1J. This is related to 
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FIG. 1: The particle yield ratios described by the present multi-component hadron gas model. The best fit for 
y/s]^ = 4.9 GeV is obtained for T ~ 131 MeV, fi B ^ 539 MeV, pia - ~ 16 MeV ( left panel), whereas for 
ysiw = 17 GeV (right panel) it is obtained for T ~ 147.6 MeV, \i B ~ 218 MeV, /i /3 ~ -2.1 MeV. A yield ratio of 
two particles is denoted by the ratio of their respective symbols. 



an essential improvement of the kaons and their ratios in the present model. Thus, the A anomaly [U [30] 
is not seen at this energy. 

In the SPS energy range we used only the NA49 mid-rapidity data for all ratios. There are two main 
reasons for such a selection. First, the NA49 are self-consistent and have relatively small error bars for 
all energies. Second, it is well known that right these data are traditionally the most difficult ones to be 
described within the thermal model (TIIT1 [8} 131} 132] . Therefore, in order to demonstrate the new possibilities 
of the multi-component hadron resonance model we concentrate on the NA49 data fitting. In contrast to 
PQ, we included into the fit procedure Q/ir~ and 3/tt~ ratios, but excluded from it the dependent H/A 
and f2/H ratios for hyperons. The results for the highest SPS energy -^/snn = 17.3 GeV are shown in the 
right panel of Fig. [TJ These results are compared to the NA49 mid-rapidity data for pions, kaons and 
(anti)protons IMIEI], for the set of strange (anti)hyperons [35H37] and for <\> meson [38] . 

The variation of the R n and Rk radii immediately allowed us to notably improve the description of 
K + /tt + ratio and all the ratios involving the strange hyperons and pions (for instance, look at E - /-7r~ 
and Q/ir~). This also led to a slight improvement of K~ /K + ratio. However, a slight change for the 
total strangeness of kaons means a larger change of the strange hyperons densities. Although the most 
problematic ratios at this energy, namely A/A, are J2/fi, are improved only marginally compared 

to [1], but as one can see from the right panel of Fig. [TJ the crossed ratios of H/A and fi/H, which were not 
fitted, are automatically reproduced well. The obtained results for the chemical FO temperature T and 
baryonic chemical potential \xb almost coincide with the values T ~ 152 MeV, [is — 226 MeV found in [1] 
for this energy for the fitting the NA49 data alone. However, the resulting quality of our fit at this energy 
of collision is essentially better: x 2 /dof ~ 1.57 determined in this work against x 2 /dof ~ 2.78 found in 
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The same trend is seen for all the SPS energies: a small variation of kaon hard-core radius and a 
vanishing pionic hard-core radius systematically improve the fit quality of the NA49 data. After such a 
fitting of the most 'hard' data it is clear that a high quality description of other data sets (or of all data 
sets) is possible, if the experimental data of different collaborations are reanalyzed and become consistent 
with each other. 

Since the RHIC high energy data of different collaborations agree with each other, we just analyzed 
the STAR results for ^/s^F = 130 GeV |39H12] and 200 GeV [Hffi]. For the main subject of the present 
work the exotic ratios are not of a great importance and, hence, for the RHIC energies we fitted the same 
set of hadronic ratios as for the highest SPS energy (see the right panel of Fig. [I]). The variation of the 
pionic and kaonic hard-core radii practically does not affects the fit quality at the RHIC energies. This is 
clearly seen from the comparison of the best fit results for y^s/vw = 130 GeV found here T ~ 163.1 MeV, 
HB ^ 27.3 MeV and the values T ~ 162.5 ± 5.5 MeV, [i B — 35 i 11 MeV found in [1J for a combined 
fit of PHENIX and STAR data assuming that the pion data do not contain any contribution from weak 
decays. The better agreement is seen for the chemical FO parameters at ^snn = 200 GeV obtained here 
T ~ 162.2 MeV, fi B ^ 17.6 MeV and the values T ~ 160.5 ± 2 MeV, \l b - 20 ± 4 MeV found in [1J for a 
combined fit of PHENIX and STAR data excluding p/ir~ and <j)/K~ ratios. 
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FIG. 2: y/sNN dependence of the chemical FO temperature (left panel) and baryonic chemical potential (right panel) 
found here within the multi-component model (circles). For a comparison the corresponding quantities for a model 
with a single hard-core radius R = 0.3 fm [5] (squares) are also shown. 

The above results are not surprising, since from Fig. [2] it is seen that the present multi-component fit 
almost reproduces the values of the chemical FO temperature and baryonic chemical potential obtained 
in the model with a single hard-core radius for all hadrons [5] which in its turn very well reproduces the 
results of [lj for hadronic yield ratios. However, the main result obtained within the present model, the 
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FIG. 3: ^snn dependences of K + /tt + (left panel) and A/ir~ (right panel) ratios obtained here within the multi- 
component model are compared to that ones found within the one component model [5]. 

best description of the Strangeness Horn, i.e. K + /n + ratio, alone with A/tt~ ratio, is demonstrated in 
Fig. [3j These two ratios were thoroughly studied within the one component hadron resonance gas model 
in [1] and then their analysis (together with other ratios) was continued in [Sj. As we mentioned earlier for 
the one component model the fit of ratios is not affected by the excluded volume correction. If, however, 
one imposes some additional constraints like the artificial baryon charge conservation criticized in [5j, then 
some small correction (below 5 %) can appear. This is the reason why the one component model fit of 
K + /n + and A/n~ ratios found in [IJ [8] for a single hard-core radius R = 0.3 fm is essentially worse than 
the one component model fit with the same hard-core radius shown in Fig. [2j Evidently, in Refs. [H [8] 
the fit of A/tt~ ratio indicates a problem with too steep rise as a function of -^/snn compared to the 
data, while both of these ratios show too slow decrease compared to the data. Note that such a behavior of 
K + /tt + and A/tt~ ratios is typical for almost all statistical models (see Figs. 7 and 8 and the corresponding 
discussion in |7j). Evidently, the too steep rise in A/tt~ behavior is a consequence of the A anomaly [T| I30j 
discussed above. The one component fit of the present approach does not indicate such difficulties for 
A/ir~ ratio, while the slow decrease in K + /tt + ratio still is there. However, the multi-component approach 
really removes such a defect in K + /ir + without spoiling the other ratios including A/tt~ one as it is seen 
from Fig. [2] The results for other ratios are available and will be shown in a longer work. 

Actually the best fit for K + /tt + ratio found here practically coincides with the dashed curve drawn in 
Fig. 4 of [H] which was obtained assuming an existence of the Hagedorn mass spectrum of hadrons |45j . 
However, we do not share this hope of the authors of Ref. [8j. Also we do not agree with such an estimate 
and, hence, we cannot accept it as a real solution of the puzzling problem. Our skepticism is based on 
the following facts. First of all, we note that inclusion of the hypothetical states with the masses up to 3 
GeV should essentially modify not only K + /ir + results, but all other ratios in uncontrollable way. Then 
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the authors of [5] would spoil their own results reported in many nice works. Second, it is unclear to us 
why in this case one should stop at hadron mass of 3 GeV and not to increase it till 10 GeV or even to 
100 GeV? It is clear that already in the former case the one component excluded volume description would 
lead to a whole complex of problems that are typical to the Hagedorn mass spectrum and in order to get 
rid of them one will have to unavoidably introduce the excluded volume which is proportional to a mass 
of a heavy resonance. Such models are well known |46} I47j . but then in the framework of such models 
the heavy hadronic resonances should be regarded as quark-gluon bags and the whole treatment should be 
completely changed. 

Furthermore, in [8] the estimates based on the Hagedorn mass spectrum inclusion were done without 
accounting for the large resonance width and without knowing the branching ratios of these hypothetical 
resonances. Although the mass dependence of width of heavy resonances was found within the finite width 
model of quark gluon bags |46tl47| and recently it was successfully verified on the Regge trajectories of heavy 
mesons |48j . the possible channels of their decays and the corresponding branching ratios are completely 
unknown yet. However, the more serious issue is that the finite width model \46\ [4"7] explains that the 
huge deficit of the empirical hadronic spectrum compared to the Hagedorn one is due to the fact that the 
'missing' hadrons with the masses above 2.5 GeV and with the large width are the quark gluon bags, which 
are extremely suppressed (by fifteen-sixteen orders of magnitude) compared to the stable hadrons up to 
the temperatures of about half of the Hagedorn temperature. Therefore, we again come to a conclusion 
that the 'missing' hadrons should not be included into the hadron resonance gas model spectrum, but they 
should be attributed to the spectrum of quark gluon bags. The practical consequence out of these facts is 
as follows: due to the short life time r ~ J ^0.5 fm/c of the bag of mass M > Mq ~ 2.5 GeV, by the 
time of chemical FO such bags should have been, probably, completely decayed into the stable hadrons 
and light hadron resonances. Since up to the moment of chemical FO the chemical equilibrium is assumed 
to exist, then all thermodynamic quantities of pions (or other particles appeared from these bags) should 
(locally) have their equilibrium values in accordance with the hadron resonance gas model spectrum, i.e. 
at chemical FO the result of the bag decays should not bee seen. 



IV. CONCLUSIONS 



In the present work we considered the hadron resonance gas model with multi-component hard-core 
radii, i.e. we treat the pion R n and kaon Rk hard-core radii as independent fitting parameters compared 
to the hard-core radius of baryons Rb = 0.5 fm and that one of all other mesons R m = 0.3 fm. Such 
an approach allows us to essentially improve the quality of the global fit of hadron yield ratios derived 
from the mid-rapidity data measured at the AGS, SPS and two highest RHIC energies. Thus, for R^ = 
fm and Rk = 0.35 fm we found x 2 /dof ~ 1.018 which is the best value for the global fit compared to 
other analyses. It is necessary to stress that at SPS energies the present fit included only the NA49 data, 
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which are usually hard to be fitted by the hadron resonance gas model. The suggested approach allows us 
to drastically improve the quality of the data description and, as a consequence, it is able to completely 
describe the Strangeness Horn irregular behavior for the first time. Thus, the developed approach gives a 
simple solution to the puzzle of the Strangeness Horn description without the need to use the hypothetical 
hadron resonances of masses up to 3 GeV which, so far, are not observed in the experiments. 

The found small values of hard-core radii are consistent with the results of analysis done in [12] that up to 
the temperatures of about 170 MeV the hadron-hadron repulsive and attractive interaction contributions 
into the system pressure practically compensate each other. In fact, we determined the second virial 
coefficients of hadrons using the statistical model with the multi-component hard-core repulsion. Evidently, 
such an approach can be used to further improve the description of other particle ratios. Thus, the suggested 
multi-component model provides us with a practical way to extract the second virial coefficients for all 
hadrons and tabulate them as a function of temperature as this is done for usual gases. The main problem, 
however, is related to the poor quality of existing experimental data. The present analysis clearly shows 
that to accurately determine the hadronic second virial coefficients and, thus, to provide the community 
with the data allowing in principle to extract the statistical measure of interaction for any pair of hadrons 
we need much better data up to -y/s/vTv — 20 GeV. We hope that the Dubna Nuclotron and the future 
colliders NICA and FAIR will successfully resolve at least the half of this task. 
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